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ABSTRACT 

A general quasilinear transport parameter for particle drift in arbitrary turbulence 
geometry is presented. The new drift coefficient is solely characterized by a nonresonant 
term and is evaluated for slab and two-dimensional turbulence geometry. The calcu- 
lations presented here demonstrate that fluctuating electric fields are a key quantity 
for understanding quasilinear particle drift in slab geometry. It is shown that parti- 
cle drift does not exist in unpolarized and purely magnetic slab fluctuations. This is 
in stark contrast to previous models, which are restricted to slab geometry and the 
field line random walk limit. The evaluation of the general transport parameter for 
two-dimensional turbulence geometry, presented here for the first time for dynamical 
magnetic turbulence, results in a drift coefficient valid for a magnetic power spectrum 
and turbulence decay rate varying arbitrarily in wavenumber. For a two-component, 
slab/two-dimensional turbulence model, numerical calculations are presented. The new 
quasilinear drift, induced by the magnetic perturbations, is compared with a standard 
drift expression related to the curvature and gradient of an unperturbed heliospheric 
background magnetic field. The considerations presented here offer a solid ground and 
natural explanation for the hitherto puzzling observation that drift models often de- 
scribe observations much better when drift effects are reduced. 

Subject headings: cosmic rays — drifts — turbulence 

1. INTRODUCTION AND MOTIVATION 

Particle transport in random magnetic fields plays a key role in space physics and astrophysics. 
The knowledge of turbulence properties is crucial for understanding (anisotropic) particle diffusion 
in a collisionless, turbulent and magnetized plasma such as the solar wind or the interstellar medium. 
Of particular interest are processes governing particle diffusion across an ambient magnetic field, 
and perpendicular diffusion and particle drifts appear to be the primary mechanisms. Transport 
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coefHcients describing these two processes are thus key quantities for several hehospheric and as- 
trophysical settings such as solar modulation of cosmic rays and diffusive shock acceleration, e.g., 
at the hehospheric shock of termination or at shocks triggered by super novae. 

It is generally accepted that the modulation of cosmic rays in the hcliosphcrc is described 
adequately by Parker's equation (Parker 1965) taking into account (anisotropic) spatial diffusion 
of cosmic rays by the diffusion tensor K. This tensor describes the interaction of cosmic rays with 
the turbulence. In a local orthogonal coordinate system with one axis along the direction of the 
mean magnetic field, it takes on its simplest form and reads 

K=\ -KA «±,2 1 . (1) 




The diagonal elements describe diffusion along (k||) and perpendicular (k_l,i, k±^2) to the mean 
magnetic field. The off-diagonal, antisymmetric entry is related to effects of curvature and 
gradient drift in a nonhomogencous hcliosphcric magnetic field (HMF). For an unmodified Parker 
spiral (Parker 1958), ka reads (Jokipii et al. 1977) 

i^A = (2) 

where v and Rl denote the speed and the Larmor radius of a particle, respectively. The derivation 
of equation (2) is based on the assumption that the HMF is undisturbed, i.e., electromagnetic 
turbulent fields are not taken into account. 

The importance of large-scale drifts of cosmic rays in the HMF was first recognized in the 
late seventies by Jokipii et al. (1977) and subsequently taken into account, via equation (2), in 
enumerable numerical studies (e.g., Potgieter & Burger 1990; Webber et al. 1990; Potgieter et al. 
1993; Burger &; Hattingh 1998). While being indispensable for a description of multi-dimensional 
large-scale modulation, it has been demonstrated that drift effects do not have to be included in 
all cases to reproduce observations (see, e.g., Reinecke et al. 1993; le Roux &; Fichtner 1997). 

It has been pointed out in several numerical studies that the effect on modulation, as described 
by equation (2), is too large if the numerical results are compared with observations (e.g., Potgieter 
et al. 1987, 1989; Potgieter & Burger 1990; Webber et al. 1990). To obtain a better agreement 
between simulations and observations, it was suggested to use a reduced amount of drift effects 
at low and intermediate particle energies (Potgieter et al. 1987, 1989; Burger 1990). In spite of 
its importance and the pressing need of rigorous theories, only one model has been developed in 
the past being able to explain a possible reduction of the drift at intermediate energies (Bieber & 
Matthaeus 1997). In brief, they obtain the expression 

3 1 + (tO)2 ^"^^ 

with CI being the relativistic particle gyrofrequency. The timescale r is associated with the decor- 
relation of a particle trajectory from an unperturbed helical orbit. For their calculations, Bieber 
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k, Matthaeus (1997) assume slab turbulence geometry and argue that nonresonant field line ran- 
dom walk (FLRW) is the major agent for the decorrelation of particle trajectories. Based on this, 
they demonstrate that drift effects can be suppressed at low and intermediate particle energies. 
Another model, developed earlier by Forman et al. (1974) on the basis of quasilinear theory (QLT) 
and FLRW in slab geometry, is formally based on the same expression for k^, equation (3). It 
has been shown, however, that this model fails in explaining suppressed drift effects, since their 
approach is valid only for cosmic ray energies greater than ~ 3 GcV at Earth (Bieber &; Matthaeus 
1997). Besides being valid for slab geometry only, both models take into account purely magnetic 
fluctuations and were developed for static turbulence only. Finally, for completeness, it should be 
mentioned that equation (3) is formally the same as for hard-sphere scattering in a magnetized 
plasma (Gleeson 1969), where r is then the scattering time. Furthermore, a similar structure was 
found for collisional particle transport parameter in a thermal equilibrium plasma (Balescu et al. 
1994) and expressions for cocfRcients being of similar forms were already considered in the mid 
1960s (cf., e.g., Toptygin 1985, and references therein). 

Upon looking into the literature, one becomes aware of the fact that a rigorous, quasilinear 

treatment of particle drift in an electromagnetic plasma wave turbulence with arbitrary geometry 
does not exist. The purpose of this paper is, therefore, to present a solid QLT treatment of particle 
drift in slab and, particularly, 2D turbulence geometry. 

The structure of the paper is as follows: Section 2 introduces the governing QLT equations of 
motions (Sec. 2.1) required for the evaluation of second-order velocity cross-correlation functions 
(Sec. 2.2). General Fokker-Planck coefficients and drift coefficients valid for arbitrary turbulence 
geometry and arbitrary plasma wave dispersion relations are presented in Sec. 2.3. Drift coefficients 
for an electromagnetic plasma wave turbulence with slab geometry are presented in section 3. 
Section 4 gives a QLT particle drift coefficient for 2D turbulence geometry, presented here for the 
first time. The new drift coefficient is valid for an arbitrary wavenumber dependence of the magnetic 
power spectrum and the turbulence decorrelation rate. Numerical calculations for a two-component 
turbulence model are presented, together with the conclusions, in section 5. 



Spatial diffusion and drift coefficients are commonly calculated from an ensemble of particle 
trajectories. For statistically homogeneous and stationary fluctuations, the so-called Taylor-Green- 
Kubo (TGK) formula (e.g., Kubo 1957) is often employed. 



2. QUASILINEAR FORMALISM 



t 




(4) 







in the limit t ^ oo. Here, Vi is the ith Cartesian component of the single particle velocity. The 
brackets < . . . > denote an ensemble average over the relevant distribution of particles. For large 
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coherence time ^, the second-order velocity correlation function < Vi{0)vj{$,) > must go to zero, 
and the integral in equation (4) approaches a constant value for t — ^ oo. 

Within the context of QLT, the transport coefficients Kij can be written as (Schlickeiser 2002, 
Eq. (12.3.25)) 

1 

i^ij = \j diiDxiX^, (5) 
-1 

where /x = uy/u is the pitch-angle of a particle having the velocity component v^^ along the ordered 
magnetic field Bq. The subscripts Xi and Xj denote guiding center coordinates in the ith and jih 
Cartesian direction, respectively, and Dx^Xj denotes Fokker-Planck coefficients of the form 

oo 

Dx,x, =^Jd^< Xi{0)Xj\O > ■ (6) 



They represent the interaction of a particle with electromagnetic fluctuations. The relation in 
equation (6) correspond to the TGK formula (4) used earlier by Bieber &; Matthaeus (1997) and 
Forman et al. (1974). In QLT, however, one has to perform an additional average with respect to 
II. 

To evaluate the QLT transport coefficients for particle drift in a large-scale magnetic field with 
superimposed electromagnetic fluctuations, the corresponding Fokker-Planck coefficients Dxy and 
Dyx have to be calculated. This requires the temporal variations of the quiding center coordinates. 
The particle equations of motion then enable one to calculate the corresponding second-order 
velocity cross-correlation functions. 



2.1. EQUATIONS OF MOTION 

According to Schlickeiser (2002, Eqs. (12.1.9d) and (12.1.9e)), the perpendicular components 
of the fluctuating force fields can be written as 

X{t) = -vcos(t){t)^/l^—^ + — [6Ey + ii-6BA (7) 

ijQ -DO V C / 

Y{t) = -vsm<|){t)^/^^^-p--^(6E,-^i-6By), (8) 

where ^ and c denote the gyrophase of the particle and the speed of light, respectively. Note that 

the Cartesian components SB,j.^y^^^ and ^E'^, y y of the fluctuating electromagnetic field are used and 
not the helical description, i.e., left- and right-hand polarized fields. For the further treatment of 
equations (7) and (8), a standard perturbation method is applied. To do so, it is convenient to 
replace in the Fourier transform of the irregular electromagnetic field the true particle orbit x(t) 
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by an unperturbed particle orbit, yielding 

/+00 „ 
d3MB(k,i)e*^W-* = (fk6B{k,t)JniW)exp 



in[ip — (pit)] + ik\\V\\t 



(9) 



and an analogous expression for 5E. The quantity J„(VF) is a Bessel function of the first kind 
and order n. The particle gyrophase for an unperturbed orbit is given by 4>{t) = (po — fit, where 
the random variable 00 denotes the initial gyrophase of the particle. Furthermore, the abbrevia- 
tion W = k±RL^/T — is introduced, where Rl = v/ft is the Larmor radius. The relativistic 
gyrofrequency is given hy Q = qBo/{^mc) with m being the mass and q the charge of the particle, 
7 is the Lorentz factor. The angle tp results from the wavenumber representation kx = k± cos ip 
and ky = k±simp. With equation (9), the equations of motion (7) and (8) can be manipulated to 
become 



Xit) 



oo „ 



d^kexp 



m [t/; — (p{t)] + ik\\v\\t 



(10) 



X < - 



J, 



n+1 



(T^)e^^+e-^^J„_i(ty) 



6Bii + -JJW) (SEy + H-6B, 
" V \ c 



CO „ 



m [ip — 0(t)] + ik\\V\\t 



(11) 



X < i- 



Jn+i{W)e''^-e~''l'Jn-iiW) 



5Bn--JniW) (6E^-n-6B, 



For the evaluation of equation (5) , it is convenient to consider now the nature of the electromagnetic 
turbulence. Here, the "plasma wave viewpoint" is employed and it is assumed that the turbulence 
consists of a superposition of N individual plasma wave modes, i.e., 

N N 

SB(k, t)=Y^ SB^{k) expi-iujt) ; 5E(k, t) = Y^ SW{k) exp{-iojjt). (12) 
j=i j=i 

Here, Wj(k) = uj^RCk) + ^^j(k) is a complex dispersion relation of wave mode j, where uj^r is the 
real frequency of the mode. The imaginary part, rj(k) < 0, represents dissipation of turbulent 
energy due to plasma wave damping. 

Restricting the considerations to transverse ((JE-' • k = 0) fluctuations and using Faraday's law, 
the turbulent electric field can easily be expressed by the corresponding magnetic counterparts, 
yielding 

(13) 



-^^[bB\^,-bBik 
Furthermore, it is convenient to use the Bessel function identities 



2n 

Jn-l{W) + Jn+l{W) = ^Jn{W) ; Jn-l{W) - Jn+l{W) = 2j;,(W), (14) 
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where the prime denotes the derivation with respect to W. With equations (13) and (14), Eqs. 
(10) and (11) can readily be rearranged to become 



oo „ 



exp 



J n=—oo 



m [^p — (f){t)] + i(A;||U|| — u;j)t 



(15) 



X Jn{W) 



k± II 



Y{t) 



N oo „ 



j=l n=— oo ' 



m [ip — (pit)] + «(A:||U|| — LOj)t 



X Jn{W) 



a^5B{ + bSBi 



where the following complex functions have been introduced: 



V k^ 



H. 



(16) 



(17) 



2.2. VELOCITY CROSS-CORRELATION FUNCTIONS 



Having determined the equations of motion, one can now proceed to calculate the second-order 
velocity cross-correlation functions < X(0)y*(^) > and < y(0)X*(^) > entering equation (6). The 
procedure for the calculation is relatively lengthy, but can be carried out with simple algebra. The 
calculations for both correlation functions are analogous, and the calculations are restricted to 
< X(0)y*(^) >. Multiplication of equation (15) with the conjugate of equation (16) leads to 



oo oo 



x(o)y*(0 



^ j n=—oo m=— oo 
Ajn{W)Jm(W) 



hEE E I d'k I d^k eMx) 

K: 

kxky _^ 
-aa 

± 

kx — * 



=^aa* • {5Bl5B\:) + 66* • {SBidBi*) 
.k±k± II " 



,-^7'*> k^j 



+i^af ■ (6Bl6Bi*) + ^a*b ■ (5Bi6B 
k± II " k± 



-^./T^J^{W)J^{W)^ 



^La* . {6Bl6^{') + b* ■ {6BUBI* 
k±_ 



kx 



k± lk± 



kx 



a ■ [SBiSBi*) + 6 • {SBlSBi*) 



2\ 



kykx ji 



k±k_i 



(18) 
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with X = i{n'ip — mip) — i{n — m)(f)o — i{k^^v^^ — a;^+mJ7)^. The bar notation used over some quantities 
indicates that they have to be evaluated for wavevector k and time ^. 

A further simphfication of (18) can be achieved only if an average with respect to the random 
variable the initial gyrophase of the particle, is applied. For this, the relation 

27r 

^ J d(j)oexp[i{n-m)(l)o] = Snm (19) 



is used, where Snm = ^or n ^ m and unity for n = m. Furthermore, the ensemble average is 
applied and it is assumed that the Fourier components at different wave vectors are uncorrelated. 
Introducing the subscripts a and /3 for Cartesian coordinates, the ensemble averages of the magnetic 
field fluctuations then read 

< 5Bi6W; >=< 6Bi{k)6B^;{k') >= 6{k - k')Pi^{k). (20) 

The uncorrelated state implies ip = ip, W = W and Uj = uj. Equation (18) can then be manipu- 
lated to become 

2 oo 

<X(0)y*(O> = E d\k2^eM-^{k\\v^\-'^*j+nm] (21) 

^0 j n=-oo-' 



^{JI{W) 



k:ckyaa*P^^^^ + klbb*Piy + k,,k±ab*Pf^y + kyk±a*bP^^^^ 



iVr^Jn{W)J^{W) 



Pjll {k^a* + kla) + k±kyb*Ply + k±k^bP^ 



-{l-,/)kyk, [j'n{W)fP^ 

and an analogous expression for the cross-correlation function < y(0)X*(^) > governing the trans- 
port parameter kyx was derived. Both are expressed by a sum of three individual terms, and each 
term is accompanied by a specific factor through which the components of the magnetic correlation 
tensor P^^(k, ^) enter the cross-correlation functions. 



2.3. FOKKER-PLANCK COEFFICIENTS 



Having determined the velocity cross-correlation for kxy in the previous section, one can 
now proceed and evaluate the corresponding Fokker-Planck coefficient Dxy- Upon substituting 
equation (21) into (6), one obtains 



2 °o r T? 

1 n=— oo -L 



Jl{W)Fi 



XY 



lJniW)j'niW)G' 



XY 



XY 



(22) 
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with the auxiliary functions 

Fj^Y = fca:fcj^aa*P|[|| + kibb*Piy + k:^k^ah*P^^y + kyk±a*bP-^ 



(23) 



Pjll {k^a* + kla) + k±kyb*Pf. + kMP-^ 



(24) 



H'^Y — (1 A'' y • 

The integration with respect to ^ leads to the complex resonance function, 

Tj + t{k\\v\\ — ujj^R + nO) 
T-] + {k\\v\\-Uj,R + nQ.f 



TZj = exp [— ^(^ii^'ji — ^j,R + nO,)^ + Fj^] = — 



(25) 



(26) 



which describes interactions of the particles with the plasma wave turbulence. The calculations for 
the Fokker-Planck coefficient Dyx are analogous to the calculations for Dxy and result in 



2 oo „ 

^0 r,.=-r^ J 



kl 



YX 



J n=—oo 

with the corresponding auxiliary functions 



F^^ = fc3.fcj^aa*P|j'|| + k'\bh*Pi^ + k^k±a*bP^^^ + kyk^ab*P^^^ 



(27) 



(28) 



^YX 



Pjll {kla* + kla) + kx.kybPl^^ + kM*Pl 



(29) 



jjj _ H ..2 



'YX 



{1 - fl^) k^kyP' 



(30) 



Equations (22) and (27), one of the main results of this paper and presented here in this general 
form for the first time, allow one to calculate QLT drift coefficients for arbitrary turbulence geom- 
etry, where the turbulence consists of transverse wave modes with dispersion relations depending 
arbitrarily on wavevector. 

Further treatment of the aTixiliary functions (23) to (25) and (28) to (30) requires a cer- 
tain representation for the correlation tensor P^g- Different representations for will alter the 
underlying mathematical and physical structure of the Fokker-Planck and, therefore, the drift co- 
efficients. Here, a representation is chosen commonly used in the literature. Following, e.g., Lerche 
&; Schlickeiser (2001), the nine components of P^g can be expressed as 



A^ik±,kii 



Sap - + ^0"^(fc±,A;||)ea/3i.y 



(31) 
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where denotes the magnetic helicity, S^p is Kronecker's delta and eap^ is the Levi-Civita tensor, 
is the wave power spectrum. With equation (31), one can now proceed to evaluate the Fokker- 
Planck coefficients (22) and (27). The calculations are very laborious and result in quite lengthy 
expressions. For the sake of overview, a simplification is introduced which concerns the complex 
functions a and b given by equation (17). Their complex nature results from the imaginary part 
of the wave mode dispersion relation, i.e., the dissipation rate Fj. It enters a and b via Faraday's 
law used to express the turbulent electric field components by their magnetic counterparts. For 
instance, consider the first term in equations (23) and (28). The quantity aa* can also be written 
as 




Analogously, one can cast the other expressions, such as bb* or ab*, into contributions including 
either (jJj^r or Tj. The simplification is to neglect the contributions given in terms of Fj, without 
the loss of insight or important information. This reduces all equations substantially, and under the 
assumption that plasma wave damping is weak, this simplification seems to be reasonable. Upon 
substituting the appropriate components of (31) into the auxiliary functions (23) to (25) and (28) 
to (30), one arrives at 

-2 [k^ky [{aRk± - bRk\\)'^ - bj^k'^] - ia^bRk'\k{aRk±_ - bRk\\)\ (33) 



where 



F^x = "p" i^^h [(.(^Rk± - bRk\\f - byi^] + ia^bRk'j_k{aRk± - bRk\\)] (34) 

A^k± 

^XY = V 1 - A*^ [k±iaRk± - bRk\\) - 2ia^bRkkxky] (35) 

^YX = ^-^V^^--J^[k'i{aRk±-bRk\\) + 2ia^bRkk3;ky] (36) 

H'xY = H^yx = Ani-^.')^-^, (37) 



W vk-^ vk'^ 

According to equations (22) and (27), one has to take the real part of the wavenumber integral. In 

view of the resonance function (26) and Eqs. (33) to (37), it is clear that further calculations involve 
both the real and imaginary part of equation (26). This is not the case for quasilinear perpendicular 
diffusion, where the real part of the resonance function (26) is required only (Stawicki 2004). 

In order to proceed, equations (33) to (37) arc substituted into the coefficients (22) and (27) 
and the relation aRk± — bRk\\ = {k\\^v\\ — cOj^R + nQ)/v is used. The coefficients Dxy and Dyx can 
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then be expressed as 




-N 

R sin ^ cos + { (39) 
+N 



with the operator 

representing resonant wave-particle interactions. The contribution N represents non-resonant be- 
havior and reads 



xjcT^' (j^-f^) kJliW) + Vl^k^UW)J'^{W)^. (41) 

A closer inspection of the operator (40) results in the finding that R sin ^ cos V' = 0, since no addi- 
tional dependences in ip are assumed. This implies vanishing resonant wave-particle interactions. 
Quasilinear particle drift is then solely characterized by the non-resonant term N, equation (41). 
Upon using the relation (5), one obtains 

1 

-1 

As expected for particle drift, kxy and Kyx arc antisymmetric. The nonresonant part, Eq. (41), 
allows to determine the drift coefficient (42) for different turbulence geometries. Motivated by 
theoretical work (Zank &: Matthaeus 1992) and observations (Matthaeus et al. 1990; Bieber at el. 
1996), slab and 2D turbulence geometries are assumed here, and each is considered in turn. 



3. DRIFT COEFFICIENT FOR SLAB GEOMETRY 

In slab turbulence geometry, the wavevectors are all either parallel or antiparallel to the back- 
ground magnetic field, and the wave power spectrum can be given by 

A^ = 5i(fe||)^. (43) 

To calculate kxy or, simultaneously, kyx for slab geometry, the Bessel functions in equation(41) 
are considered in the limit W oc k± ^ 0. The first term in the braces contributes only for n = 0, 
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since Jn(0) = 1 for n = only. Concerning the second term, it can be shown that it vanishes 
completely for — and all n. The pitch-angle integration is elementary and (42) then yields 




(44) 



For further process, the wave power spectrum -S's(A;||) = C{u)Xs{SBl)^{l + /cyAg) is employed. 

Here, is the bend-over scale, (SBl)'^ is the slab variance and 2u = 5/3 is the inertial range 
spectral index. Furthermore, C{u) = {2-\/n)~^r{i') /T(u — 1/2), where r(x) denotes the Gamma 
function. Concerning the plasma wave dispersion relation, it is assumed that forward (j = +1) and 
backward {j = —1) propagating shear Alfven waves with real frequency coj^r = jvAk\\ form the slab 
turbulence. For va/v « 1, the second term in equation (44) then becomes zero. This does not 
imply that the dynamical behavior of the turbulence due to dissipation is not taken into account, it 
is rather suppressed since the argument of the logarithmic expression becomes unity. Upon using a 
constant magnetic helicity, the wavenumber integration can be performed analytically and equation 
(44) yields 



KxY — ~'^YX — ~'^0 



" [Jb^J \6BsJ 



(45) 



Here, kq = TT^XsVAi^B / Bq)'^ , where < ^ = 5B1/5B'^ < 1 measures the fraction of the slab 
contribution to the total turbulent magnetic energy, 5B^ = 5Bg + SB^j;)- The evaluation of Eq. 
(42) for 2D turbulence geometry is presented in the next section. A closer inspection of the slab drift 
coefficient (45) results in the following findings: First, K^y entirely determined by the magnetic 
helicities a'^ and (t~ of forward and backward propagating wave fields, respectively. Second, 
depends on neither the charge nor the mass of the particle and is, therefore, independent of particle 
properties. 

In view of the independence of n^y of particle properties, it is instructive to recall a standard 
zeroth-ordcr drift, the E x B drift. A charged particle moves with the drift velocity v/j = c(E x 
B)/B^ if an electric force acts normal to the background magnetic field. This drift is identical for 
all charged particles and, therefore, independent of particle charge, mass and velocity. Equation 
(45) reveals the same feature. An enlighting approach for a comparison is to replace E by the 
fluctuating field (5E. The drift velocity v^j is then a random quantity and a corresponding velocity 
cross-correlation function and, therefore, drift coefficient can be derived. Based on this, it can be 
shown that n^y is indeed a result of the (5E x Bq drift, where the perpendicular force results from 
the electric component of the turbulence. 

To obtain some more insight into the drift coefficient (45) , it is convenient to recall the definition 
of the magnetic helicity It is usually defined as 

, _ {5Btf - {5Blf 
-{5Bt? + {6Bir ^ ^ 
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where 5Bl and SBr denote left-handed (LHP) and right-handed polarized (RHP) field components 
of the fluctuations, respectively (see, e.g., Schlickeiscr 2002). Now a variety of wave fields with 
different polarization states and propagation directions can be considered: 

(1) LHP wave in forward direction {a'^ = +1, SBg = 0): equation (45) yields = — '^o- A. 
RHP wave in backward direction {a^ = —l,SBg = 0) leads to the same result, i.e., a LHP 
wave in forward direction can be replaced by a RHP wave propagating in backward direction. 
If both are present, the drift becomes two times stronger, = —2ko. 

(2) RHP wave in forward direction ((t+ = —1, 5B'^ = 0): equation (45) yields K^y = +'^0- A 
LHP wave in backward direction {a~ = +1, 5B^ = 0) leads to the same result. As in case 
(1), if both types of waves are present, the drift becomes two times stronger, i.e., ^XY ~ +2^0 

(3) For equal polarization states, i.e., = a~ , one finds K^y = —i^oHcCr^ with 

^ {5B^f - {5B-f 

" {5B+y + {5B-y ^ ' 

being the normalized cross helicity (Schlickeiscr 2002). It represents the ratio of the intensities 
of forward (j = +1) to backward (j = —1) propagating wave fields and is sometimes also 
referred to as "Alfvenicity" . Obviously, K^y becomes zero for a vanishing net polarization. 
Furthermore, it changes sign if predominantly LHP or RHP wave fields are present, i.e., with 
the reversal of polarity. 

At a glance, the drift coefficient for slab geometry, Eq. (45), is solely determined by the magnetic 
helicity and the real frequency ojj^R. If one of them is neglected, no QLT particle drift occurs 
in slab turbulence. However, Forman et al. (1974) and Bieber & Matthaeus (1997) argued that 
FLRW governs drift in a static and purely magnetic slab turbulence. In stark contrast to this are 
equations (44) and (45). They clearly show that particle drift in slab geometry requires turbulent 
electric field components and is due to the net polarization only, i.e., the magnetic helicity aK 
Furthermore, the nonresonant FLRW limit used by Forman et al. (1974) and Bieber &; Matthaeus 
(1997) is based on the concept of a magnetic power spectrum at zero wavenumber, SB^{k\\ = 0). 
This limit can be achieved only for static {Tj = 0) turbulence conditions. The real part of the 
resonance function (26) yields then a Dirac delta distribution 5{k\\) required for the FLRW limit. 
However, the real part represents resonant interactions. As it is shown in Sec. 2.3, the real part is 
not important for QLT particle drift, and so FLRW. QLT drift in slab geometry is a consequence 
of the nonresonant term, Eq. (41). The latter results from the imaginary part of equation (26) and 
can not explain the FLRW limit. 



4. DRIFT COEFFICIENT FOR 2D GEOMETRY 

Having derived the drift coefficient for a plasma wave turbulence and gained the insight that 
fluctuating electric fields are required for QLT drift in slab geometry, the evaluation of the non- 
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resonant contribution (41) for 2D geometry is presented in this section. For 2D geometry, the 
wavevectors are perpendicular to the mean magnetic field, and the wave power spectrum can be 
given by 

A^{k^,^) = Si^{k^)^-^. (48) 

Obviously, shear Alfven waves as used in the previous section for slab turbulence can not con- 
tribute, since ujj^r = in equation (41) for the power spectrum (48). The calculations are, there- 
fore, restricted to purely dynamical magnetic fluctuations. Since the concept of a superposition of 
individual wave modes does not apply anymore, the j-nomenclature is dropped. The only modifi- 
cation concerns the resonance function (26). To take into account the dynamical behavior of purely 
dynamical magnetic fluctuations, Bieber et al. (1994) defined two models: the damping as well as 
the random sweeping model. For the damping model, they suggested a dynamical behavior of the 
turbulent energy being of the form 

< SB^{k,t)5B;{k,t + >=Pa/3(k)exp(-z/ee), (49) 

where fc represents decay of turbulent magnetic energy. The resonance function (26) has to be 
adapted in this respect, and the plasma wave dissipation rate Tj is simply replaced by the decor- 
relation rate z^c- 

Upon substituting equation (48) into (41), replacing Tj by Vc and applying the //-average, one 
obtains 

4". - -^E / / ".^^MWMW). (50) 

"-^ 

Here, the fact was used that the first term in equation (41) is an odd function in jj, and, therefore, 
vanishes, due to the jU-average. As a consequence, the magnetic helicity a (not assumed to be zero) 
does not influence quasilinear drift in 2D turbulence geometry. The /x-integration in equation (50) 
can be carried out analytically and the detailed calculations are presented in Appendix A. There, 
it is shown that equation (50) can be manipulated to become 

oo 

4y = -4x = -^^^ / dk±S2D{k±)I{C, z) (51) 



with 

7r/2 

7(C,^) = j d9 



sin(26l) sinh(26'z) 
^3 sinh(2;7r) cos^ 9 





2C cos e cos(2C cos 6") + (4C^ cos^ 6" - 1) sin(2C cos 6) 



(52) 



where the abbreviations C, = k±RL and z = Uc/Cl are introduced. Furthermore, Rl = v/Cl is 
the Larmor radius. Equation (51) is valid for a power spectrum and a decorrelation rate varying 
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arbitrarily in wavenumber k±. It is noteworthy that a similar expression can also be derived 
for quasilinear perpendicular diffusion in 2D turbulence geometry (Stawicki 2004). The integral 
representation (52) results from the /[/-integration and has to be evaluated for further progress. 
Unfortunately, an analytical solution to this integral does not exist, and any progress requires 
a numerical treatment. Figure 1 shows numerical computations of equation (52) as functions of 
C = k±RL for three different values of z. For illustrative purposes, z is assumed to be a constant in 
k±. Using small argument approximations for the hyperbolic sine, it can be shown that Eq. (52) 
is independent of z and a function of ( only. Furthermore, note that /((", z) and, therefore, K^y 
changes sign with the reversal of a positive to a negative particle charge state. 



2 Q _ I I I I I 




1 2 3 4 5 6 

Fig. 1. — Plot of numerical solutions to Eq. (52) representing 7(^, z) as a function of ^ = k±RL for 
three different values of z (see legend). 



The limit ^ <^ 1, implying Rl ^ kj_ , leads to an instructive, analytical solution for equa- 
tion (52). To show this, small argument approximations for the circular functions are used, i.e., 
cos(2<^cos6') ~ 1 and sin(2(^cos0) ~ 2(^cos0 and inserted into I{(,z). Partial integration then re- 
sults in /(^ ^ l,z) = 8/3(1-1-2;^). Consequently, one obtains for the drift coefficient the expression 

oo 


in the limit Rl ^ kj^, where the relation Uc = has been used. Since the correlation time 
Tc is still undetermined, equation (53) can be considered as being valid regardless of whether 
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is smaller than, larger than, or of order unity. However, in view of the restriction k±RL <S 1, it 
becomes obvious that (53) is vaUd only for low/intermediate particle energies if parameters are 
assumed being typical for the hcliosphcrc. An eyecatching feature of equation (53) is the term 
including the dimensionless product TcO. It is formally the same as those given in equation (3) 
used by Forman et al. (1974) and Bieber & Matthaeus (1997) for the FLRW limit. 

5. NUMERICAL CALCULATIONS AND CONCLUSIONS 

To demonstrate the potential and flexibility provided with the new drift coefficients, the re- 
maining wavenumber integration in Kxy^ ^1- (^1)' solved numerically and a two-component, 
slab/2D turbulence is considered. Since the individual contributions arc simply additive, the total 
drift coefficient Kp = i^xY ~^ '^XY' induced by the fluctuations, is introduced, where equations (45) 
and (51) are used for K^y '^Wy^ respectively. For the 2D component, the magnetic power 
spectrum 

S2D{kA.) = (1 - i)C{u)\2DSB''{l + klxlj,)-- (54) 

is employed. Here, as for the slab drift coefficient, C(i/) = (2y^)~^r(z/)/r(i/ — 1/2) and < ^ = 
dBg/{6Bg +5i?|^) < 1, where is the total variance of the 2D-component. The corresponding 

bend-over scale is given by X2D and 2i^ = 5/3 is the inertial range spectral index. For the numerical 

treatment, the turbulence dccorrclation rate Vc has to be specified entering equation (51) via z = 
Vf-jVl. In analogy to Bieber et al. (1994), is assumed to be of the form = avAk± where the 
parameter < a < 1 allows adjustment of the strength of dynamical effects. The case a = 
represents the magnetostatic limit, a = 1 describes a strongly dynamical magnetic turbulence. 

For the numerical computations, conventional parameters being typical for the heliosphere are 
used. The ratio {SB/Bq)^ is assumed to be 0.2, until otherwise noted. The parameter a is set to 
unity and the Alfven speed va is chosen to be 50 km s"^. The background magnetic field Bq is 
given by 4 • 10"^ Gauss. For all calculations, it is assumed that the turbulent magnetic energy has 
only a small fraction in its slab component (say 20%) and is dominated by the 2D turbulent energy 
(80%), yielding = 0.2. Until noted otherwise, it is assumed that = IOA2D = 0.03 AU (Bieber 
et al. 1994). For convenience, all Figures show absolute values of the drift coefficients. 

Figure 2 shows numerical result of k'xy^ Eq- (51), as a function of the proton Larmor radius 
Rl normalized to the slab bend-over scale As for ^ = 0.2 and three different values of \2d/^s (see 
legend). The ratio Rl/^s is proportional to the particle rigidity R. For the computations, it is 
assumed that the slab contribution is unpolarized, i.e., (T+ = (7~ = 0, yielding K^y = 0- The 
additional solid line visualizes the large-scale drift coefficient ka, Eq. (2). A closer inspection of 
Figure 2 results in the following finding: and ka obey the same power law dependence in 

Rl/Xs or, alternatively, rigidity R for R^/Xg <C 1, namely R^. This is expected. As it is shown 
in Sec. 4, the general 2D drift coefficient (51) reduces to equation (53) for Rl <^ kj^^. Obviously, 
KA and then reveal the same variation in Rl/Xs or, synonymously, R. However, independent 
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of the ratio Md/^s-, i^xy becomes a constant for Rl/Xs ^ 1, whereas ka scales with R. For 
completeness, Figure 3 illustrates numerical results of for protons, electrons and helium for 
X2d/\s = 10. 




10-3 10-2 10-1 100 101 102 

Fig. 2. — Numerical solutions to Eq. (51) representing function of Rl/\s for protons 

and three different values of the ratio A2D/A5 (see legend). The slab component is suppressed by 
assuming an unpolarized state ((7+ = a~ = 0), and the ratio SB"^ /Bq is chosen to be 0.2. It is 
assumed that (5-B|^ : 5B1 = 8:2. The additional solid line visualizes the standard large-scale drift 
coefficient ka, Eq. (2). 

As explained in Sec. 1, is valid for an unperturbed, unmodified Parker spiral only (see 
Eq. (2) and the comments following it), while Kp describes effects due to the two-component 
turbulence. Generally, it is expected that additional (electro) magnetic turbulent fields alter drift 
effects, particularly for low and intermediate particle energies. Whether with or without turbulence, 
particle motion is affected by curvature and gradient drift effects. The standard coefficient ka might 
therefore be considered as the limit of a more general drift coefficient for 5B — ^ for which 
vanishes. Since the fluctuating fields are superimposed to the heliospheric background magnetic 
field, it is assumed that individual drift effects induced by Bq and 5B are simply additive (at least in 
the local orthogonal coordinate system). This results in a total drift coefficient kt = ka — Kp- Note 
that both KA as well as change signs with the reversal of the particle charge state. According to 
Eq. (51), but also equation (45) for the slab contribution, one has to take into account an additional 
minus sign, resulting in the difference of ka and Kp- The slab contribution is independent of the 
particle charge state and depends solely on the polarization of slab turbulence. In this respect, the 



/ 



Protons 



/ 



/ 



Electrons 



/ 



10-3 - 



/ 



/ 



Helium 



10-4 




J I 



J I 



J I 







10-3 



10-2 



10-1 



100 



10» 



102 




Fig. 3. — Plot of numerical solutions to Eq. (51) for 5B^ /Bq = 0.8, = IOA2D and three different 
particle species: protons, electrons and helium. 

magnetic hclicity of the slab component entering equation (45) and the normalized cross liclicity 
He, Eq. (47), are quite uncertain parameters. A rigorous theory describing wavenumber and radial 
variations of the slab helicities does not exist. Usually, for the 2D component, He is assumed to be 
zero for heliocentric distances beyond ~ 1 AU (e.g., Zank et al. 1996), but this might probably not 
be the case (see Matthaeus et al. 2004). 

Figure 4 shows results for (solid line), ka (dotted line) and their difference, kt (dashed 

line), for the turbulence level 5B^/Bq = 0.8 and an unpolarized slab contribution, i.e., /tf-y = 0. 
Here, protons are considered. At a glance, the drift coefficient ka is substantially reduced at low and 
intermediate particle rigidities. At very low rigidities, kt is almost two magnitudes smaller than 
the standard drift ka, but reveals the same power law behavior in rigidity, i.e., R^. With further 
increase in R, kt varies as Ry^, rolling over to Kp oc R at high rigidities, where the large-scale 
relation ka dominates. 

As mentioned in Sec. 1, a reduction of the amount of drift effects at low and intermediate 
rigidities was suggested earlier by, e.g., Potgieter et al. (1987) and Burger (1990) and was then con- 
sidered theoretically by Bieber & Matthaeus (1997). Assuming SB'^/Bq ~ 1, Bieber Sz Matthaeus 
(1997) found a i?^ power law behavior for the scaling of drift effects at intermediate energies. How- 
ever, absolutely central to their approach is the assumption that FLRW governs the reduction of 
the drift. This was not assumed here. The FLRW limit is rather excluded, since the imaginary, 
nonresonant part of the function (26) governs quasilinear particle drift, and not the real part. 
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Fig. 4. — Numerical solution of (solid line) for protons and SB^/Bq = 0.8 as well as = 
IOA2D. The slab component is suppressed, n^y ~ assuming an unpolarized state, i.e., 

= a~ = 0. The dashed curve represents the total drift coefficient kt = ka — hf- 

Furthermore, the argumentation by Bieber & Matthaeus (1997) is valid for slab geometry only. 
For the numerical computations given in Figure 4, the slab component is explicitly excluded by 
assuming an unpolarized state. Obviously, the reduced amount of drift effects shown in Figure 4 
results rather from the presence of the 2D component than from slab turbulence and, least at all, 
FLRW. Depending on the ratio 5B'^/Bq, the influence of magnetic perturbations on drift effects, 
induced in the local orthogonal coordinate system by curvatures and gradients of the global and 
unperturbed background magnetic field, leads to a significant reduction of these drift effects. The 
result BJ/'^ for SB'^/B^ = 0.8 is relatively close to the result by Bieber & Matthaeus (1997), i.e., 
R^. However, Figure 4 shows that kt recovers to continue with R^ with decreasing rigidity or ratio 
Rl/Xs- This is different when slab turbulence is polarized. 

The influence of the slab contribution is shown in Figure 5. Here, the same linestyle and 
parameters are used as those for Figure 4, but a net polarization of the slab component is assumed. 
For this, a relatively weak right-hand polarization of forward and backward propagating wave fields 
is chosen, i.e., = a~ = —0.1 (see Sec. 3: item (3), page 12). For the variances of the wave fields, 
the ratio {SB')"^ : (6B~^)'^ = 9 : 11 is assumed, implying for the normalized cross helicity Hc = 0.1. 
This implies a slightly larger amount of turbulent energy in forward than in backward propagating 
wave fields. A comparison of Figures 5 and 4 results in the insight that the slab drift coefficient 
^XY becomes important only at very low values of R^/Xg, i.e., very low particle rigidities. For this 
range, n^y exceeds by an order of magnitude and kt becomes a constant for a slab component 
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Fig. 5. — Same as Fig. 4, but the slab component is assumed to be weakly right-hand polarized 
= a~ = —0.1). A slightly larger amount of turbulent energy in forward propagating wave 
fields is assumed, i.e. He = 0.1. The dashed-dotted curve represents — K^y, Eq. (45). The dashed 
curve represents the total drift coefficient kt = ka — kf- 



being weakly right-hand polarized. 

The reduction in the amount of drift effects and the change in rigidity dependence varies 
for different particle species. Figure 6 shows numerical results similar to those given in Figure 
5, but now for electrons. First, as expected, the 2D drift coefficient varies as on R at low 
and intermediate values of Rl/Xs, indicating the relativistic nature of the electrons at such low 
rigidities. As a consequence of this, kt scales with R^/"^ instead of i?^/^ and rolls then over to 
kt (X R at high rigidities. 



The aforementioned long-standing assumption of a vanishing normalized cross helicity he for 
heliocentric distances beyond ~ 1 AU (e.g., Zank et al. 1996) was made for the 2D turbulence com- 
ponent only. Unfortunately, a similar treatment for the slab component does not exist. Assuming 
that the same holds for slab turbulence, it implies that vanishes in the outer heliosphere and 
is present only within Earth's orbit. In this case, the 2D component dominates QLT particle drift 
throughout the outer heliosphere, since it does not depend on any helicity (see Eq. (50) and the 
comments following it). However, using a more advanced transport model including cross helicity, 
Matthaeus et al. (2004) recently relaxed the assumption of a vanishing cross helicity for 2D turbu- 
lence. Their study indicates that cross helicity might even be present at around 10 AU. Assuming 
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Fig. 6. — Numerical solution of (solid line) for the same parameters and linestyle as those 
used for the computations shown in Fig. 5, but here for electrons. 



that the same statement holds for slab turbulence, this would imply that k^y might become im- 
portant not only at Earth, but also for heliocentric distances within Jupiter's and Saturn's orbit, 
especially for low rigidity particles and provided that slab turbulence is polarized (even if weakly 
only). 

Another very interesting issue, which can be addressed in the context of magnetic helicity and 
its importance for QLT particle drift in slab geometry, concerns the helicity density of the large-scale 
HMF. Bicber et al. (1987) have considered the topological properties of the Parker interplanetary 
field and have shown that the helicity density of the Parker field is negative north of the heliospheric 
current sheet and positive south of the current sheet, independent of the sign of the solar poloidal 
magnetic field. Bieber et al. (1987) argue that the magnetic helicity of interplanetary small-scale 
turbulence may well be related to the helicity of the large-scale Parker field. This implies here 
that the magnetic helicity and, therefore, the polarization state of the slab component north of 
the heliospheric current sheet is opposite in sign to the polarization of the slab component south 
of the current sheet. In view of this implication, the slab drift coefficient K^y, Eq. (45), would 
reverse sign across the heliospheric current sheet, regardless of the particle charge state. Since 
the 2D drift coefficient equation (51), is independent of the magnetic helicity, only the drift 

coefficient for slab geometry would be affected by the change of sign across the current sheet. 
However, changes sign with the reversal of the particle charge state (see Sec. 4). The above 
conclusions would add new elements to possible influences of the magnetic helicity (polarization) on 
heliospheric cosmic ray transport and their solar modulation, but the investigation of their impact 
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on cosmic ray solar modulation is far beyond the scope of this paper. 

Finally, the question arises if the heliospheric transport of the so-called pick-up ions might not 
be affected by, at least, a polarized slab turbulence. Effects resulting from curvature and gradient 
drifts and spatial diffusion are usually assumed to be small and, therefore, negligible for this low 
energy particle population (see Rucihski et al. 1993). 

I thank R. A. Burger, J. Minnie, H. Moraal and M. S. Potgieter for valuable and inspiring 
discussions. Support by the South African National Research Foundation (NRF) is acknowledged. 

A. DERIVATION OF THE 2D DRIFT COEFFICIENT 

To derive the drift coefficient Eq. (51), for 2D turbulence geometry, the identity 



uwmw) = ^^^^y^^ [jum - ji^,{w)] (ai) 



4n 

is used. Equation (50) can then be cast into the form 



2^ _ -kvRl 



°° ^ oo 



^ J dk^SMk±) j c//x(1 - /x^) -^r^2 [^n-i(W^) - Jl^xm\ , (A2) 

^ „ n— 1 







where z = vd^. Employing the relations (Gradshteyn & Ryzhik 1965) 

7r/2 



Jl^^{W) = ^(-1)"^^ j dOJoilW COS e)cos{2e[n± I]), (A3) 





one obtains 

•7r/2 



Jl-i{W) - Jl+i{W) = ^(-1)""^ j deJoi2Wcos9) sin(20n) sin(2e). (A4) 



Upon substituting equation (A4) into (A2) and making use of formula (1.445.4) of Gradshteyn & 
Ryzhik (1965), i.e., 

y-/ .y»-i »sin(2gn) _ tt sinh(2gz) 

^ ^2 + ^2 2 sinh(z7r) ' ^ ' 

n=l ^ ' 

equation (A2) can be manipulated to become 

oo 7r/2 

W = -^i?L / (ifcx^^T^ / (iesin(2^)sinh(2ez)/^(^,C) (A6) 
J smh(z7r) J 



22 



with 



1 

i^{dX) = j d^i{l-^i^)Jo{2Ccose^/l-fi^). 



(A7) 



Here, W = kxRi^/T^ = CV^ 

— iJ? is used. The //-integration can be solved analytically 
(Gradshteyn & Ryzhik 1965) to obtain 



Ji/2(2Ccos^) ^ J3/2(2Ccose) 



2 [ (2Ccos0)V2 (2Ccos^)3/2 _ 
2C cos e cos(2C cos e) + {AC cos^ ^ - 1) sin(2C cos 6) 
8(3 cos3 e '' 



(A8) 



where Ji_,_^(x) denotes spherical Bessel functions of the first kind. The drift coefficient (A6) can 
then be expressed as 



-^Rl J dk^S2Dik±)I{(:,z) 



(A9) 



with the function z) as defined in equation (52). 
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